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Origin of the material dependence of T c in the single-layered cuprates 
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In order to understand the material dependence of T c within the single-layered cuprates, we study 
a two-orbital model that considers both d x 2_ y 2 and d z i orbitals. We reveal that a hybridization 
of d z 2 on the Fermi surface substantially affects T c in the cuprates, where the energy difference 
AE between the d x 2_ y2 and d z 2 orbitals is identified to be the key parameter that governs both 
the hybridization and the shape of the Fermi surface. A smaller AE tends to suppress T c through 
a larger hybridization, whose effect supersedes the effect of diamond-shaped (better-nested) Fermi 
surface. The mechanism of the suppression of d-wave superconductivity due to d z 2 orbital mixture is 
clarified from the viewpoint of the ingredients involved in the Eliashberg equation, i.e., the Green's 
functions and the form of the pairing interaction described in the orbital representation. The 
conclusion remains qualitatively the same if we take a three-orbital model that incorporates Cu 4s 
orbital explicitly, where the 4s orbital is shown to have an important effect of making the Fermi 
surface rounded. We have then identified the origin of the material and lattice-structure dependence 
of AE, which is shown to be determined by the energy difference AEa between the two Cu3d orbitals 
(primarily governed by the apical oxygen height), and the energy difference AE V between the in- 
plane and apical oxygens (primarily governed by the interlayer separation d). 

PACS numbers: 74.20.-z, 74.62.Bf, 74.72.-h 
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I. INTRODUCTION 

Despite the fact that the history of the high-T c 
cuprates exceeds two decades, there remain a number 
of fundamental questions which are yet to be resolved. 
Among them is the significant variation of T c among 
various materials within the cuprate family. It is well 
known that T c varies strongly with the number of Cu0 2 
layers, but an even more basic problem is the T c vari- 
ation within the single-layered materials. This is high- 
lighted by La 2 _ a; (Sr/Ba) a ;Cu04 with a T c ~ 40 K ver- 
sus HgBa 2 Cu04+,5 with a T c ~ 90K, with a more than 
factor of two difference despite similar crystal structures 
between them. 

Empirically, it has been recognized that the mate- 
rials with T c ~ 100K tend to have "round" Fermi 
surfaces, while the Fermi surface of the La system is 
closer to a diamond shape, and this has posed a long- 
standing, big puzzle, since, the latter would imply a rel- 
atively better nesting [l], 0. The materials with rounded 
Fermi surfaces conventionally have been analyzed with a 
single-band model with large second (i 2 (> 0)) and third 
(t^(< 0)) neighbor hopping integrals, while the "low-T c " 
La system has been considered to have smaller t 2 ,t3. 
This, however, has brought about a contradiction be- 
tween theories and experiments. Namely, while some 
phenomenological[l] and t-J model [3, @ studies give a 
tendency consistent with the experiments, a number of 
many-body approaches for the Hubbard-type models with 
realistic values of on-site interaction U show suppression 
of superconductivity for large t 2 > and/or £3 < 0, as 



we shall indeed confirm below[6[. 

To resolve this discrepancy, we have introduced in 
ref.p} a two-orbital model that explicitly incorporates 
the d z 2 orbital as well, while the usual wisdom was that 
the d x 2_ y 2 orbital suffices. The former component has in 
fact a significant contribution to the Fermi surface in the 
La system. We have shown that the key parameter that 
determines T c is the hybridization of the two-orbitals, 
which is in turn governed by the level offset AE between 
the d x 2_ y 2 and d z 2 Wannier orbitals. Namely, the weaker 
the d z 2 contribution to the Fermi surface, the better for 
d-wave superconductivity, where a weaker contribution 
of the d z 2 results in a rounded Fermi surface (which in 
itself is not desirable for superconductivity), but it is the 
"single-orbital nature" that favors a higher T c supersed- 
ing the effect of the Fermi surface shape. Recently, there 
have also been some other theoretical studies regarding 
the role of the d z 2 orbital played in the cuprates [8h11||. 

Purpose of the present paper is two-fold: By elabo- 
rating the two-orbital model, we investigate (i) why the 
d z 2 hybridization on the Fermi surface suppresses the su- 
perconductivity, and (ii) what are the key components 
that determine the material dependence in the level off- 
set between d x 2_ y 2 and d z 2. In examining point (ii), in 
addition to La2CuC>4 and HgBa2Cu04+a considered in 
ref.Q, we also construct effective models of the single- 
layered cuprates Bi 2 Sr 2 Cu06 and Tl 2 Ba2Cu06 to reveal 
how these materials can be classified in terms of the cor- 
relation between the lattice structure parameters and the 
level offsets of various orbits. 
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II. CONSTRUCTION OF THE TWO-ORBITAL 
MODEL 

A. Band calculation 

Let us start with the first-principles band 
calculation [12] of La 2 Cu04 and HgBa 2 Cu04, whose 
band structures are displayed in FigUJ The lattice 
parameters adopted here are experimentally determined 
ones for the doped materials [l3l [Hj]. In both cases, 
there is only one band intersecting the Fermi level. 
Therefore, the d x 2_ y 2 single-orbital Hubbard model, 
or the Cu-d^-y? 



0-po- three-orbital model (whose 



antibonding band crosses the Fermi level) has been 
adopted in conventional theoretical studies. A large 
difference between the two materials in the shape of the 
Fermi surface is confirmed in Figj2] As mentioned in 
the Introduction, the materials with a rounded Fermi 
surface have been modeled by a single-orbital model 
with large second f£2(> 0)] and third [ts(< 0)] neighbor 
hopping integrals [lj. It has been noticed that when the 
fluctuation exchange approximation (FLEX) [lH fl6| is 
applied to this model, a rounder Fermi surface coming 
from larger second and third neighbor hoppings results 
in a suppressed T c , as we have shown in Fig.l of ref|7|. 
A calculation with the dynamical cluster approximation 
(DCA) shows that a negative £2 works destructively 
against d-wave superconductivity [l7| . and a more re- 
alistic DCA calculation that considers the oxygen p a 
orbitals for the La and Hg cuprates also indicates a 
similar tendency fl8j. 



La 2 Cu0 4 



HgBa 2 Cu0 4 




FIG. 1: First-principles band structures of La2Cu04 (left) 
and HgBa2Cu04 (right). The top (bottom) panels depict the 
strength of the d x 2_ y 2 (d z 2) characters with the radius of the 
circles. 



B. The two-orbital model 

To resolve the above problem for the d x 2_ y 2 single- 
orbital model, we now focus on other orbital degrees of 
freedom. In fact, FigfT] shows that in the La system the 
main band has a strong d z 2 character around the N point 
on the Fermi surface that corresponds to the wave vectors 
(ir, 0), (0, 7r) in a square lattice. This has been recognized 
from an early stage of the study on the cuprates fl ffi22l . 
and more recently, it has been discussed in refs.[l|, |23[ 
that the mixture of d z 2 character to the main component 
determines the shape of the Fermi surface. Namely, the 
large d z 2 contribution in the La system makes the Fermi 
surface closer to a square (i.e., a diamond), while in the 
Hg cuprate the d z 2 contribution is small and the Fermi 
surface is more rounded (as confirmed in the following). 

In order to understand the experimentally observed 
correlation between the Fermi surface shape and T c , 
we consider a two-orbital model that takes into ac- 
count not only the d x 2_ y 2 Wannier orbital but also 
the d z 2 Wannier orbital explicitly 0j- A first-principles 
calculation [l^, H3 is used to construct maximally local- 
ized Wannier orbitals [1^, [26j], from which the hopping 
integrals and the on-site energies of the two-orbital tight- 
binding model for the La and Hg cuprates are deduced. 
Thus obtained band structures of the two-orbital model 
for the La and Hg cuprates are shown in Figf31 along 
with the Fermi surface for the band filling of n = 2.85 
(n =number of electrons per site) , which corresponds to 
0.15 holes per Cu atom. 

In the present two-orbital model, the d x 2_ y 2 Wan- 
nier orbital originates primarily from the Cu 3d x 2_ y 2 
and the in-plane O 2p a orbitals. On the other hand, 
the d z 2 Wannier orbital originates mainly from the Cu 
3d z 2 and apical O 2p z orbitals. Namely, this model in- 
corporates two types of d-p a antibonding states, where 
the former spreads over the Cu02-plane while the latter 
along the c-axis(Fig0]). Table|Tj shows the parameter val- 
ues of the present model, from which we can identify that 
d x 2_ y 2 — d z 2 interorbital hopping occurs mainly between 
nearest-neighbor Cu sites, which gives rise to the orbital 
mixture. Because the d x 2_ y 2 — d z 2 hopping integrals are 
similar for the La and Hg compounds, the onsite energy 
difference AE = E x 2_ y 2 — E z 2 between the two orbitals 
can be used as a measure of the d z 2 mixture. Note that 



La2Cu04 



HgBa 2 Cu0 4 





FIG. 2: (Color online) The Fermi surface of the La2Cu04 
(left) and HgBa2Cu04 (right) with 0.15 holes/Cu atom. 
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the interorbital hoppings have different signs between x 
and y directions, i.e., the matrix element has the form 
— 2ti[cos,(k x ) —cos(k y )], so that the d x 2_ y 2 — d z 2 mixture 
is strong around the wave vectors (n, 0), (0,7r) (N point 
in the La cuprate), while small around \k x \ = \k y \. 

In Tabl^Jwe also show the parameters for the single- 
orbital model obtained by the similar method. In the 
single-orbital model, the "d x 2_ y i" Wannier orbital effec- 
tively contains the d z 2 orbital in the tail parts of the 
Wannier orbital. 



C. Correlation between the curvature of the Fermi 
surface and AE 

The d z 2 orbital contribution has also a large effect on 
the curvature of the Fermi surface [ll, [23[, which can in- 
deed be seen from Tabled] as follows. In the single-orbital 
model, the La cuprate has smaller t% and t% as compared 
to the Hg cuprate (with the ratio ( | i 2 1 + |t3|)/|ii| being 
0.14(0.37) for La(Hg)), resulting in the smaller curvature 
of the Fermi surface in the former as mentioned. On the 




HgBa 2 CuQ 4 





\ 












m 



































r m x r z 



FIG. 3: (Color online) The band structure (with Ef = 0) in 
the two-orbital (d x 2_ y 2-d z 2) model for La2Cu04 (left column) 
and HgBa2Cu04 (right). The top (middle) panels depict the 
weights of the d x 2_ y 2 (d z 2) characters with thickened lines, 
while the bottom panels are the Fermi surface for the band 
filling of n — 2.85. The inset shows the band structure of 
the three-orbital model (see text) for La system, where the 4s 
character is indicated by thick lines. 



TABLE I: Hopping integrals within the d x 2_ y 2 orbital for 
the single- and two-orbital models (upper half), interorbital 
hopping (middle), and AE = E x 2_ y 2 — E z 2 (bottom). 





1-orbital 




2-orbital 






La 


Hg 


La 


Hg 


t(d x 2_ y 2 -S> d x 2_ y 2) 

ii[eV] 


-0.444 


-0.453 


-0.471 


-0.456 


t 2 [eV] 


0.0284 


0.0874 


0.0932 


0.0993 


t 3 [eV] 


-0.0357 


-0.0825 


-0.0734 


-0.0897 


(M + M)/I*i| 


0.14 


0.37 


0.35 


0.41 


t{d x 1_y1 ^ d z l) 

h[eV] 






0.178 


0.105 


fe[eV] 






small 


small 


t 3 [eV] 






0.0258 


0.0149 


AE[eV] 






0.91 


2.19 
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FIG. 4: (Color online) The top panel shows the main compo- 
nents of the two Wannier orbitals (having different types of 
a bonding) considered in the present two-orbital model. The 
bottom panel shows the definition of the level offsets AE, 
AE d , and AE V . 



other hand, in the two-orbital model that considers the 
d z 2 orbital explicitly, the ratio (1^1 + I^D/I^il within the 
d x 2_ y 2 orbital changes to 0.35(0.41) for the La(Hg). The 
value is nearly the same between the single- and two- 
orbital modelling of Hg, while the value is significantly 
increased in the two-orbital model for La. The reason 
why £2 and £3 in the two-orbital model for La is large 
as compared to that in the single-orbital model can be 
understood from Figf5]as follows. Let us consider the di- 
agonal hopping (f 2 ). There is a direct (d x 2_ y 2 — d x 2_ y 2) 
diagonal hopping, but there is also an indirect diago- 
nal hopping that becomes effective when AE is small, 
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that is, d r 



In the single-orbital 



model, where the d z 2 component is effectively included 
in the d x 2_ y 2 Wannier orbital, the contribution of the 
d x 2_ y 2 — > d z 2 — > d x 2_ y 2 path is effectively included in 
ti . The latter contribution has a sign opposite to that of 
the direct diagonal hopping (the reason of which will be 
clarified later) , so that we end up with a small effective £2 
in the single-orbital model when AE is small as in the La 
cuprate. A similar argument applies to £3. Conversely, 
the Hg cuprate has a large AE so that the d z 2 contri- 
bution barely exists in the single-orbital model, and the 
ratio (|i2| + I^D/I^il is similar to that in the two-orbital 
model. 

In the La cuprate, the d x 2_ y 2 and d z 2 orbitals strongly 
mix around the N point, so that the upper and lower 
bands repel with each other there, and the saddle point 
of the upper band that corresponds to the van Hove sin- 
gularity is pushed up to nearly touch the Fermi level for 
the band filling of n — 2.85. Thus the Fermi surface al- 
most touches the wave vectors (n,0), (0, 7r). In the Hg 
cuprate, there is no such splitting of the two bands, and 
the saddle point stays well below the Fermi level, result- 
ing in a rounded Fermi surface that is closed around the 
wave vector (jr, n). 



(a) Single-orbital Model 




(b) Two-orbital Model 



(c) Three-orbital Model 



FIG. 5: (Color online) Origin of the effective second-neighbor 
hopping (t-2 in the single-band model, (a)) in the two- 
orbital(b) and three-orbital(c) models. 



III. MANY-BODY CALCULATION OF THE 
SUPERCONDUCTIVITY 

A. Calculation method 

We now consider a many-body Hamiltonian based 
on the two-orbital tight-binding model discussed above, 
which is given, in the standard notation, as 



H = EEE f 



E E E c Lct c j 



E ( U E n *Mt n */4 + U ' E E n iP° n W"' 



J 



9 E E c h 



H<r C iH<T' ' c iva' c ivo 



(1) 



where i,j denote the sites and /i, v the two-orbitals, 
while the electron-electron interactions comprise the in- 
traorbital repulsion U, interorbital repulsion U', and 
the Hund's coupling J(= pair- hopping interaction J' ). 
Here we take U = 3.0 eV, U' = 2.4 eV, and J=0.3 
eV[13|. These values conform to a widely accepted, first- 
principles estimations for the cuprates that the U is 7-10t 
(with t ~ 0.45 eV), while J, J' ~ 0.1U. Here we also ob- 
serve the orbital SU(2) requirement U' = U — 2J. 

To study the superconductivity in this multi-orbital 
Hubbard model, we apply FLEX approximationfl^, [l6|, 
[28j]. In FLEX, we start with the Dyson's equation to 
obtain the renormalized Green's function, which, in the 
multi-orbital case, is a matrix in the orbital representa- 
tion as G/j^/j, where l\ and I2 are orbital indices. The 
bubble and ladder diagrams consisting of the renormal- 
ized Green's function are then summed to obtain the spin 
and charge susceptibilities, 



Xs(q) 



X°(q) 



1-Sx°(q)' 



Xc(q) = 



X°(q) 



1 + Cx°(q) 

where q = (q,iu> n ), the irreducible susceptibility is 
X? u i 2 ,i 3 ,u(l) = ^2 G hh( k + q) G uU k )> 



(2) 



(3) 



(4) 



with the interaction matrices 



£1/2 ,^3^4 



h = h = h = I4 
h = h ^ h = h 
h = h 7^ h = h 
h — I4 7^ h = 



(5) 
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C; 



U l\ — I2 — ^3 — I4 

-U' + J W = h^h = h 

2U' -J, l 1 = l 2 ^l 3 = U 

J' l\ = Z4 7^ I2 = h- 



(0) 



With these susceptibilities, the fluctuation-mediated ef- 
fective interactions are obtained, which are used to calcu- 
late the self-energy. Then the renormalized Green's func- 
tions are determined self-consistently from the Dyson's 
equation. The obtained Green's functions and the sus- 
ceptibilities are used to obtain the spin-singlet pairing 
interaction in the form 

v s (q) - \suo)S -\cuq)c + \(S + C), (7) 

and this is plugged into the linearized Eliashberg equa- 
tion, 



AA„,(fc) = E *W(9) 

X G hh (k - q)A hu (k ~ q)G hh (q - k). (8) 



The superconducting transition temperature, T c , corre- 
sponds to the temperature at which the eigenvalue A 
of the Eliashberg equation reaches unity, so that A at 
a fixed temperature can be used as a measure for T c . 
In the present calculation, the temperature is fixed at 
k^T = 0.01 eV, which amounts to about 100K, and 
the band filling (number of electrons/site) is set to be 
n = 2.85, which corresponds to 0.85 electrons per site 
in the main band, namely, around the optimum doping 
concentration. We take 32 x 32 x 4 /c-point meshes and 
1024 Matsubara frequencies. 



why we have T c positively correlated with AE simultane- 
ously with the roundness of the Fermi surface that is also 
positively correlated with AE. This should lead to the 
experimentally observed correlation between the shape of 
the Fermi surface and T c [l], 0. 



C. Effects of the inter-orbital electron-electron 
interaction 

Thus the next important question is: why does the 
mixture of the d z 2 orbital on the Fermi surface suppress 
superconductivity? To investigate the origin, we have 
varied the interaction values to examine the strength of 
the spin fluctuations and the superconducting instabil- 
ity. The strength of the spin fluctuation is measured by 
the antiferromagnetic Stoner factor, which, for a multi- 
band system, corresponds to the largest eigenvalue of the 
matrix Sxo- 

In the result in TablCD we can compare the cases for 
U' = eV and U' = 2.4 eV, which shows that the 
strength of the spin fluctuation becomes smaller when 
U' is turned off. This should be because U' hinders four 
electrons (two d z 2 and two d x i_ y i) to come on the same 
site. Despite this, it can be seen that A is not much 
affected by U' , probably because the suppression of su- 
perconductivity due to the increased charge/orbital fluc- 
tuations (which is unfavorable for singlet d-wave pairing) 
and the enhancement due to the increased spin fluctu- 
ations roughly cancel with each other. We have also 
examined how the Hund's coupling J affects supercon- 
ductivity. A comparison between J = and J = 0.3 
shows that superconductivity is slightly suppressed when 
we turn on J, which is consistent with an observation that 
the Hund's coupling tends to suppress spin-singlet pair- 



B. Correlation between T c and AE 

Let us now investigate how the d z ? orbital affects su- 
perconductivity by hypothetically varying AE from its 
original value 0.91 eV (shown in Tabk-[I| to 4.0 eV for the 
La cuprate to single out the effect of AE. The eigenvalue 
of the Eliashberg equation A calculated as a function 
of AE in FigJS] shows that A initially increases rapidly 
upon increasing AE, then saturates for AE > 3eV. This 
means that the mixture of the d z i orbital on the Fermi 
surface around the wave vectors (jr, 0), (0, n) does in- 
deed strongly suppress superconductivity in the original 
La system, while for large enough AE the system es- 
sentially reduces to a single-orbital model, where d z i or- 
bital no longer affects superconductivity. As mentioned 
above, the d z i orbital mixture makes the Fermi surface 
more square shaped, which in itself favors superconduc- 
tivity as mentioned in sectior llL^I (e.g., Fig.l of ref0). 
Thus we can see that the effect of the d z i orbital mix- 
ture supersedes the effect of Fermi surface shape, and 
T c is primarily determined by the former. This explains 
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FIG. 6: (Color online) The eigenvalue, A, of the Eliash- 
berg equation for d-wave superconductivity plotted against 
AE = E x 2_ y 2 — E z 2 for the two-orbital (red open circles) 
and three-orbital (red filled circles) models for La2Cu04. 
Corresponding eigenvalues for HgBa2Cu04,Bi2Sr2Cu06 and 
Tl2Ba 2 Cu06 are also indicated. 
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ing. Nevertheless, the effect of J is overall small. The 
conclusion here is that the effect of the interorbital inter- 
actions on superconductivity is small, so that the main 
origin of the suppression of superconductivity is the mix- 
ture of the d z 2 orbital on the Fermi surface, which will 
be elaborated in the next subsection. 

TABLE II: FLEX result for the eigenvalue of the Eliashberg 
equation A, and the Stoner factor for various values of the 
interorbital interactions U' and J, for fixed U = 3.0 eV and 
J' = 0.30 eV 



U'[eV) 


J[eV] 


Stoner 


A 


2.4 


0.3 


0.979 


0.279 


2.4 


0.0 


0.978 


0.335 


0.0 


0.3 


0.925 


0.291 


0.0 


0.0 


0.958 


0.309 



D. Origin of the suppression of superconductivity 
by the d z 2 mixing 

Here we pinpoint why the d z i orbital component mix- 
ture degrades ei-wave superconductivity. In FigJTl we 
show the squared orbital diagonal and off-diagonal el- 
ements of the Green's function matrix spanned by the 
orbital indices at the lowest Matsubara frequency. We 
compare them for two cases; the original La cuprate and 
a hypothetical case where we increase AE to the value 
for Hg, where the hopping integrals are tuned to retain 
the shape of the Fermi surface to that of the La cuprate. 
In the hypothetical case, the interaction values are re- 
duced (U = 2.1 eV, J = J 1 = 0-117 and U' = U - 2 J) 
so as to make the maximum value of the pairing inter- 
action in the d x 2_ y 2 channel (Vim) to be roughly the 
same as that in the original La case. Then the eigen- 
values of the Eliashberg equation at T = 0.01 differ as 
much as A = 0.28 and 0.88 for the original La and the 
hypothetical cases, respectively. Let us analyze the ori- 
gin of this difference. Here we denote the d x 2_ y 2 and d z 2 
orbitals as orbitals 1 and 2, respectively. In the origi- 
nal La, compared to the hypothetical case, (i) the d x 2_ y 2 
diagonal element |Gn| 2 is smaller especially around the 
wave vectors (n,0)/(0,n), (ii) the d z 2 diagonal element 
|G 22 | 2 is much larger, and (hi) there is a substantial off- 
diagonal element |Ci2 1 2 due to the strong d z 2 orbital mix- 
ture. If we turn to the pairing interaction matrix, again 
at the lowest Matsubara frequency, in FigJSl the diago- 
nal elements have similar maximum values between the 
two cases because the interaction is reduced in the hy- 
pothetical one as mentioned above. In the original La, 
the off-diagonal element of the pairing interaction V1221 
is large compared to the hypothetical case, and the in- 
teraction is broadly peaked around (0,0). On the other 
hand, the d z 2 diagonal interaction V2222 is finite but has 
a small momentum dependence. Considering the above, 
the dominant contributions to the Eliashberg equation 
regarding the d x 2_ y 2 orbital component of the gap func- 



tion An is roughly given as 
AA U (A)~ 

- Vun(Q)Gu(k - Q)A n (k - Q)G n (Q - k) 

- Vi22i(0,0)Gai(fc)Aii(*c)G 2 i(-fc) 

- ^2222(<7)G 22 (fc - <7)A 22 (fc - q)G 22 (q - k), (9) 

where Q — (tt, 7t). If we consider a wave vector k near 
(7T,0) on the Fermi surface that has a positive An(fc), 
An (A; — Q) will be negative for the d-wave gap. Then the 
first term on the right-hand side will be positive but small 
in the original La compared to the hypothetical case be- 
cause of the small Gu especially around (n, 0)/(0, tt). 
This is the main reason why A is reduced in the origi- 
nal La compared to the hypothetical case. In addition, 
the second term, which cannot be neglected when the 
d z 2 mixture is significant, actually has a negative sign, 
and also acts to suppress A and hence T c . The inter- 
action V 2222 has small momentum dependence, so that 
this term has small contribution for a d-wave gap when 
summed over q. 

In the above comparison, we have reduced the interac- 
tions in the hypothetical case so as to make the maximum 
pairing interaction Vim nearly the same as in the orig- 
inal La. The reason we fix the strength of the pairing 
interaction is because the maximum value of the pairing 
interaction actually does not differ very much upon in- 
creasing AE in the results given in FigO The reason for 
this, despite the Fermi surface nesting becoming worse 
as we increase AE, is mainly two-fold: (i) the d z 2 orbital 
mixture on the Fermi surface becomes weaker, and (ii) 
inclusion of the self-energy in the FLEX weakens the role 
of the Fermi surface nesting played in the development 
of the spin fluctuations. Regarding the second point, in 
the random phase approximation where the self-energy 
is not considered, the Fermi surface nesting effect on the 
strength of the spin fluctuations, hence the pairing inter- 
action, is so strong that A does not increase with AE as 
in FigJfJ] (and thus the T c difference between La and Hg 
cuprates discussed later cannot be explained), although 
the effect of the increase in Gn due to the reduction of 
the d z 2 mixture is present. This may be regarded as con- 
sistent with a recent result obtained with the functional 
renormalization grou p, where the self-energy correction 
is not considered fill. 



E. d x 2_ y 2 + d z 2 + s three-orbital model 

So far we have analyzed the two-orbital model that 
considers the d„ 2_„2 and d z 2 Wannier orbitals. Actu- 
ally in ref . [ij, [23J , it has been pointed out that the " ax- 
ial state" that contains not only Cu d z 2 and ap i ca rp 2 
orbital but also the Cu-4s orbital is important in deter- 
mining the shape of the Fermi surface. In the present 
two-orbital model, the Cu 4s orbital is effectively incor- 
porated in both of the d x 2_ y 2 and d z 2 Wannier orbitals. 
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FIG. 7: (Color online) Contour plots and side views of the 
diagonal and off-diagonal elements of the squared Green's 
function for the original La and the hypothetical cases. The 
subscripts 1 and 2 stand for the d x 2_ y 2 and d z 2 orbitals, re- 
spectively. 
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FIG. 8: (Color online) Diagonal and off-diagonal elements of 
the pairing interaction depicted against (q x ,q y ). 



Namely, the Wannier orbitals have Cu 4s components 
in their tails. In order to examine the effect of Cu 4s 
orbital more explicitly, let us consider in this section a 
d x 2_ y 2 + d z 2 + s three-orbital model which takes into ac- 
count the Cu-4s Wannier orbital on an equal footing. 

In this model, the 4s Wannier orbital is a mixture 
mainly of Cu 4s and O p a orbitals. The O p a orbitals 
contain not only the in-plane Opi an ej?o- but also the api- 
cal ap i ca ip 2 . The main band originating from the 4s 
orbital for the La system is shown in the inset of Fig[3] 
While the 4s band lies well(~ 7 eV) above the Fermi 
level, the 4s orbital still gives an important contribution 
to the Fermi surface shape. Here again we estimate the 



ratio (|t2| + |<3|)/|ti| within the d x 2_ y 2, where we find a 
much smaller value of 0.10 against 0.35 in the two-orbital 
model. This means that the large ti and £3 within the 
d x 2_ y 2 Wannier orbital in the two-orbital model is mainly 
due to the d x 2_ y 2 — > 4s — > d x 2_ y 2 hopping path (Fig[5j 
bottom panel), as pointed out in ref.[l|. Then, from the 
viewpoint of the three-orbital model, £2 and £3 in the 
single-orbital model of La cuprate are small because the 
d x 2_ y 2 — } 4s — } d x 2_ y 2 and d x 2_ y 2 — y d z 2 — y d x 2_ y 2 
contributions nearly cancel with each other. The two ef- 
fective hoppings have opposite signs because the d z 2 level 
lies below d x 2_ y i while 4s lies above. 

Now we apply FLEX to this three-orbital model, where 
we vary AE = E x 2_ y 2 — E z 2 and calculate the eigenvalue 
of the Eliashberg equation as we did in subsectior illl Bl 
Here, we fix the on-site energy difference E s — Ed 2 at 
its original value when we vary AE, because the three- 
orbital model for the Hg compound has roughly the same 
E s — Ed 2 as that of the La compound. 

The result is displayed in Fig. [5] as marked with " 3- 
orbital" . We recognize that in the small AE regime the 
eigenvalue A rapidly increases with AE as in the two- 
orbital model. In the large AE regime, however, A tends 
to decrease rather than to saturate. In this regime, the 
4s level comes too close to the Fermi level, and strongly 
deforms the Fermi surface. Nonetheless, considering that 
even in the case of the Hg compound with a larger AE 
as will be discussed later, AE (3-orbital model) is still 
~ 2 eV, i.e., such a suppression of superconductivity due 
to the 4s level coming too close to the Fermi level is not 
expected in real materials. 

Thus we can conclude on the 4s orbital that, while this 
orbital has an important effect on the shape of the Fermi 
surface, the effect can be included in the two-orbital 
model, so that the FLEX results for the two- and three- 
orbital models are similar as far as the T c — AE relation is 
concerned (unless we consider unrealistically large AE). 
This is natural in that the level offset E x 2_ y 2 — E z 2 is 
smaller (~ 1 eV) than the electron-electron interaction(~ 
3 eV), while the E s — E x 2_ y 2 is much larger (~ 7 eV). 
Hence the 4s orbital can effectively be integrated out be- 
fore the many-body analysis, while the d z 2 orbital can- 
not. In this sense the two-orbital(rf 2 .2_ i) 2 — d z 2) model 
suffices for discussing the material dependence of the T c 
in the cuprates. 



IV. MATERIAL DEPENDENCE OF AE 

We have seen that the mixture of the d z 2 com- 
ponent strongly affects superconductivity, making T c 
positively correlated with AE. To further endorse 
this, we have plotted in Fig [HI the eigenvalue A 
for the two-orbital models for single-layered cuprates 
Bi 2 Sr 2 Cu0 6 [2|, Tl 2 Ba2CuO 6 [30|, and HgBa 2 Cu0 4 as 
well, whose lattice structures are shown in FigO We can 
see that these materials also fall upon reasonably well 
on the correlation between A and AE. Thus the next 
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fundamental question in understanding the material de- 
pendence of T c is: which key factors determine AE. This 
section precisely addresses that question. 




So let us first focus on how the apical oxygen height 
ho affects AEd- Namely, we construct a model that con- 
siders all of the Cu 3d and O 2p orbitals(five 3d + 3 x 4 
2p = 17 orbitals) explicitly, exploiting maximally local- 
ized Wannier orbitals, and then estimate the on-site en- 
ergy difference between Cu d x 2_ y 2 and Cu d z ? orbitals 
as AEd- We note that this AEd is something different 
from AE defined for the effective two-orbital model we 
have considered, since we now explicitly consider the oxy- 
gen 2p orbitals. In FigfTUl we plot AEd as a function of 
ho, where we hypothetically vary the height for the La 
system from its original value 2.41 A to 2.90 A . The re- 
sult shows that AEd and ho are linearly correlated. We 
have also constructed similar d-p models for the Bi, Tl 
and Hg systems, and we can see that the AEd values for 
these materials, also included in the figure, roughly fall 
upon the linear correlation for the hypothetical La sys- 
tem, which indicates that AEd is primarily determined 
by ho ■ Such a correlation has also been found in a recent 
quantum chemical calculation (Toj. where the d x 2_ y 2-d z 2 
level splitting evaluated there corresponds more closely 
to the present AEd rather than AE. 

Having seen that ho governs AEd, we next look at AE 
and the eigenvalue A in the two-orbital models for the La 
cuprate with hypothetically varied ho- As expected, AE 
in Fig llll (b) monotonically increases with ho- Then A 
(Fig[TT](a)) increases with ho, which is in accord with the 
positive correlation between AE and A discussed above 
(FigfTT] (c)). Thus ho is shown to be one of the key 
parameters that determine AE and thus T c . 

However, if we plot the corresponding values for the 
Hg cuprate, also displayed in the figure, we find that AE, 
and thus A, are larger than those for the hypothetical La 
cuprate for the same apical oxygen height between the 
two cuprates. This implies that ho and AEd are not the 
sole parameters that determine AE and hence T c , and 
another factor should be lurking. 



FIG. 9: (Color online) Lattice structures of La2Cu04, 
BiaSraCuOe, Tl 2 Ba 2 Cu06, and HgBa 2 Cu0 4 . 



A. Crystal-field effect 

Since the main components of the Wannier orbitals 
in the two-orbital model are the Cu 3d x 2_ y 2 and Cu 
3d z 2 orbitals, the crystal-field splitting between these or- 
bitals, denoted as AEd here, should be the first key fac- 
tor governing AE. Namely, materials with a larger apical 
oxygen height above the Cu02 plane(/io) should have a 
larger crystal- field splitting, [20] so that AEd, and thus 
AE, should be larger (FigfJ}. Indeed, the La compound 
has smaller ho = 2.41 A and AE, while the Hg com- 
pound has larger ho = 2.78 A and AE. 
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FIG. 10: (Color online) AEd plotted against ho- Solid (red) 
circles connected by a line is the result for the hypothetical 
lattice structure of La cuprate, while values for Bi, Tl and Hg 
cuprates are also shown. 
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FIG. 11: (Color online) The eigenvalue of the Eliashberg 
equation A (circles) when ho is varied (a) or AE(ho) is varied 
(c) hypothetically in the lattice structure of La cuprate. Also 
plotted is AE(ho) against ho(b). Diamonds in green indicate 
the values for HgBa2CuC>4. 



that of La, namely, the 0-2p z energy level with respect 
to the in-plane 0-2p CT level is much lower in Hg. 

The difference mainly comes from the crystal structure 
where the apical oxygen in the La cuprate is surrounded 
by other apical oxygens belonging to the neighboring lay- 
ers, while in Hg those oxygen atoms are much further 
apart as seen in Figj9] This gives a clue to understand the 
reason why the hypothetical La cuprate with the same 
ho as Hg has smaller AE and A; although AEd is similar 
between the two systems, AE p very much differs. Thus 
the difference between La and Hg can be attributed to 
the distance between neighboring Cu02 layers that is af- 
fected by the lattice structure, i.e., body-centered tetrag- 
onal (bet) vs. simple tetragonal. However, a similar vari- 
ance in the layer distance can occur even within similar 
lattice structures. La, Bi, and Tl compounds all have the 
bet structures, so naively one might expect similar values 
of AVa- However, AV^'s for Bi and Tl are much larger 
than that for La. This is because in Bi (Tl) there is a 
Bi-0 (Tl-O) layer inserted between the adjacent Cu02 
layers (see FigFJJ, resulting in a large Cu02 layer sepa- 
ration. 

So let us focus on the separation between the neigh- 
boring CuC>2 planes, which will be denoted as d here. 
Figure CC^a) plots AV^° c) against d for La, Hg, Tl, 

and Bi cuprates. Here we have defined AV^"^ as the 
contribution to AVa coming from the apical oxygens. 
These Madelung potentials are calculated by placing 
point charges at atomic positions, as was done in ref . [3l| . 
We have also plotted the total AVa for the four mate- 
rials, which indicates that AVa is roughly governed by 
AV^° C \ which in turn is mainly determined by d. We 
also plot AE p against d in FigfT^b) for the four cuprates. 
Here again, AE p is obtained using the model that con- 
siders the Cu-3d and 0-2p orbitals explicitly. From these 
we can see that both AVa and AE p are primarily corre- 
lated positively with the layer separation d. This in turn 
implies that AE p and AVa in FigfLWc) are positively 
correlated as well. 



B. Oxygen-orbital effects 

The above observation has motivated us to look more 
closely into the effects of oxygen orbitals. As shown in 
FigHl the Wannier orbitals in our two-orbital model the 
Cu-3dj.2 _ y 2 and 3d z 2 orbitals strongly hybridize with the 
in-plane O 2p a and apical oxygen O 2p z orbitals, respec- 
tively. Thus we can surmise that AE should also be af- 
fected by the energy difference (denoted as AE p ) between 
the in-plain p a and the apical oxygen p z . By definition, 
one can expect that AE p is positively correlated with 
AVa, the Madelung potential difference between O p i anc 
and O a pi ca i introduced by Ohta et al. as an important 
parameter that controls the material dependence of the 
T c 31]. In fact, AVa for Hg is about 7 eV larger than 



C. Classification of materials by AEd and AE P 



We have seen that AEd and AE p are mainly deter- 
mined by ho and d, respectively. Combining these, we 
can summarize the dependence of AE on the material 
and lattice structure as 



AE~f(AE d (h ),AE p (d)), 



(10) 



where / is a certain function. For instance, La and Bi 
have smaller AEd reflecting smaller ho , while Hg and Tl 
have larger AEd due to larger ho- Namely, the latter 
group tends to have larger AE. On the other hand, Bi, 
Tl, and Hg have larger d than La, so that they have larger 
AVa- We can summarize all these into a classification of 
materials in terms of AEd and AE p as a numerical ta- 
bl dllli and a kind of "phase diagram" in Fig|T3] Apart 
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FIG. 12: (Color online) (a) AV { f c) (circles) and AVa (dia- 
monds) plotted against the layer separation d for La, Bi, Hg, 
and Tl cuprates. (b) The level offset, AE P , between the in- 
plain p a and the apical oxygen p z against the layer separation 
d. (c) The correlation between AVa and AE P . 



from the effect of ho (or AEd), AE is positively corre- 
lated with AE p and thus with AVa , so that AVa and T c 
should be roughly correlated. In this sense, the so-called 
Maekawa's plot (Fig. 2 of ref.|31]) is consistent with the 
present Figj6l Also, a negative correlation between the 
occupancy of holes with p z — d z i character and T c has 
been found in ref . [22| , which is again consistent with the 
present view. 



TABLE III: The values of AE d (along with ho), AE P (along 
with d), and AE for La, Bi, Hg and Tl cuprates. 
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FIG. 13: (Color online) AE plotted against AE P and AE d 
for the four single-layered cuprates considered here. An 
oblique plane indicates a rough correlation between AE and 
(AE p ,AE d ). 



V. DISCUSSIONS 

A. Validity of the present model 

In the present study, we have adopted the LDA to de- 
rive the kinetic-energy part of the model Hamiltonian. 
The LDA calculation neglects some of the electron cor- 
relation effects, and our standpoint in the present study 
is that the remaining part of the electron correlation is 
dealt with in the FLEX calculation. One might sus- 
pect, however, that there might remain electron corre- 
lation effects, that are not taken into account in the 
present approach but can affect the accuracy of the eval- 
uation of the level offset AE between d x 2_ y 2 and d z 2 
Wannier orbitals. Our view on this point is the follow- 
ing. First, it is an experimental fact that the La cuprate 
has a square-like Fermi surface, while the Bi cuprate a 
rounded one [J- This is accurately reproduced in the 
LDA, which strongly suggests that the d z 2 component 
is indeed strongly mixed around (n, 0), (0, n), i.e., AE is 
small, in the La cuprate. Secondly, a detailed quantita- 
tive difference in AE will not affect the present conclu- 
sions. To see this, we have performed an LDA+U calcula- 
tion to obtain the kinetic-energy part of the Hamiltonian, 
varying U from to 6 eV. For La, the considerable d z 2- 
character around (tt, 0), (0, 7r) persists even at U = 6, and 
the band that intersects the Fermi level is only slightly 
changed, although AE somewhat increases with U. On 
the other hand, for Hg AE is greatly enhanced by U, but 
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this does not significantly affect the d x 2_ y 2 main band, 
since the d z 2 character is already absent at U = 0. Apply- 
ing FLEX to these LDA+U models will result in a double 
counting of the electron correlation effects because FLEX 
takes account of the first-order terms, but if we took, for 
the sake of comparison, the obtained AE, we would find 
that a considerable difference in A between La (A ~ 0.5) 
and Hg (~ 0.8) is still present even if we adopt the mod- 
ified values of AE. 



B. Possibility of higher- T c materials 

A consequence of our study is that superconductiv- 
ity in the single-layered cuprates is optimized when the 
system has a single-band nature. In such a case (as in 
Hg cuprate), the Fermi surface is rounded due to the ef- 
fect of the Cu 4s orbital. As mentioned in sectior llTAl 
the square shaped Fermi surface would be more favor- 
able for superconductivity for single-orbital systems. For 
this very reason even the Hgl^CuCi+a is not fully op- 
timized as a single-layered material. Indeed, the hypo- 
thetical La cuprate having a large AE but with a Fermi 
surface similar to that in the original La gives a larger 
A in the Eliashberg equation as we have seen in section 
IIIIDI So we have a bit of a dilemma, since it would be 
difficult to get rid of the effect of the Cu 4s orbital as far 
as the cuprates are concerned. Conversely, however, we 
can seek for other materials in which the 4s orbital is not 
effective. An example is a single-band system consist- 
ing of d xy orbitals, where the hybridization between d xy 
and 4s orbitals is forbidden by symmetry. In fact, a pos- 
sible way of realizing a single-band d xy system has been 
proposed in ref. 32]. Provided that such a system has the 
band width and the electron-electron interaction strength 
similar to those in the cuprates (since too strong or too 
weak a correlation will degrade superconductivity), it can 
possibly give even higher T c . 

VI. CONCLUSION 

To summarize, we have studied a two-orbital model 
that considers both d x 2_ y 2 and d z 2 Wannier orbitals in 



order to pinpoint the key factors governing the mate- 
rial dependence of T c within the single-layered cuprates. 
We conclude that the d z 2 orbital mixture on the Fermi 
surface is significantly degrades superconductivity. Since 
the energy difference AE between the d x 2_ y2 and the 
d z 2 governs the mixture as well as the shape of the Fermi 
surface, we identify AE as the key parameter in the ma- 
terial dependence of T c in the cuprates. Since the mixing 
effect supersede the effect of the Fermi surface nesting, a 
small AE results in a suppression of T c despite a square 
shaped Fermi surface. AE is then shown to be deter- 
mined by the energy difference AEd between the two 
Cu3d orbitals, and the energy difference AE p between 
the Opi anc p a and O ap i ca i p z , both of which are affected 
by the lattice structure. AEd is a crystal field splitting, 
which is mainly determined by the apical oxygen height, 
while AE p is found to be primarily governed by the inter- 
layer separation d. The materials that have highest T c s 
within the single-layered cuprates, Hg are Tl systems, 
indeed have AE large enough to make them essentially 
single-band. On the other hand, there is still room for 
improvement if we can suppress the effect of the Cu 4s 
mixing that makes the Fermi surface rounded, which may 
be realized in non-cuprate materials with U similar to the 
cuprates in magnitude. 
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